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Abstract 

We consider the collisions of large nuclei using the theory of McLerran and Venu- 
gopalan. The two nuclei are ultra-relativistic and sources of non-abelian Weizsacker- 
Williams fields. These sources are in the end averaged over all color orientations 
locally with a Gaussian weight. We show that there is a solution of the equations of 
motion for the two nucleus scattering problem where the fields are time and rapidity 
independent before the collision. After the collision the solution depends on proper 
time, but is independent of rapidity. We show how to extract the produced gluons 
from the classical evolution of the fields. 
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1 Introduction 


Nucleus-nucleus collisions at ultra-relativistic energy have long been recognized as 
an environment where hot dense matter is formed |l|-]3|. It has been conjectured that 
in such an environment one might produce and experimentally study a quark-gluon 
plasma Theoretical studies of quark-gluon plasma formation have typically 
assumed some initial conditions at some time after the collision was initiated, and 
then evolved the matter distributions forwards in time according to the equations 
of perfect fluid hydrodynamics |3|,|^. 

While such an approach may work well for the late stages of the collision when 
the particles are not so energetic, it does not work well for the earliest stages of the 
collision. In the earliest stages, the quarks and gluons emerge from their quantum 
mechanical wavefunction and cannot be described as a perfect fluid until at least 
enough time has passed for there to be scattering. 

In the earliest stages of the collision, the quark and gluon interactions should be 
most energetic. Such scatterings are therefore more easy to experimentally probe 
as they presumably induce hard experimental signatures which are more easily dis¬ 
entangled from backgrounds due to soft final state processes. During the hydrody¬ 
namic expansion, typically the scale of energy in the interaction is softer and more 
difficult to disentangle from backgrounds. 

There has been recent progress in attempting to describe the early evolution of 
matter produced in nuclear collisions Q . In the parton cascade model of Geiger and 
Muller, one takes the experimentally measured distribution functions for quarks and 
gluons and assumes that they may be treated as an incoherent beam of particles 
arising from each nucleus. The scattering of partons from partons is computed 
making reasonable assumptions about quantum coherence and time dilation effects. 
The system is thereby evolved from very early times in the collision until a later 
time when hydrodynamics may be applicable. In such a theory, the hard scattering 
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signals are computed and may be compared with experiment. 

The parton cascade model while elegant and well motivated, in our opinion still 
lacks some theoretical underpinning. In particular, the issue of quantum coherence 
in the initial state is treated phenomenologically, and needs a deeper understanding. 
This problem has at least two important aspects. 

The first and most glaring problem is that the partons arise from a quantum 
mechanical state. In such a state the uncertainty in momentum, Ap, times the 
uncertainty in position, Ax, is close to saturated, 

ApAx ~ 1 (1) 

For example, in the longitudinal momentum distribution of partons, the wee partons 
have a longitudinal momentum of order of lOO’s of MeV. This corresponds to a 
longitudinal size of order of fractions of a Fermi. On the other hand, in the parton 
cascade one assumes knowledge of both the position and momentum of the partons, 
since the partons are described by classical phase space distribution functions. While 
this should be true later in the collision as the scale of spatial gradients becomes 
larger, early in the collision it is most certainly violated. 

Although in the parton cascade, the assumptions on the initial distributions 
are plausible, they can at best give a qualitative agreement with precise results 
which include the effects of coherence, and at worst totally ignore some classes of 
interference phenomena. For example, one obvious problem is that for a single 
nucleus, the partons will spread out since they are an incoherent distribution of 
partons with different momentum. After some time, one therefore no longer has a 
spatially compact nucleus. 

Another class of phenomena which is not fully treated in the parton cascade 
model is the problem of coherent addition of the color charges of quarks and gluons. 
Such coherent addition is for example responsible for Debye screening, and presum¬ 
ably magnetic screening, which will serve as a cutoff for divergent transport cross 
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sections in parton-parton processes. In the parton cascade, a low momentum cutoff 
is introduced by hand, and of course results for many processes depend upon this 
cutoff. 

While the parton cascade may lack precision in many detailed computations, 
it nevertheless is outstanding for its qualitative predictions. We nevertheless would 
like to put this model onto firmer foundations, and understand clearly its limits of 
applicability. 

To begin to tackle this problem, one must understand at least some aspects 
of the quantum mechanical wavefunction of the quarks and gluons in the nuclear 
wavefunction. In the past, one rarely considered the nuclear wavefunction, and the 
structure functions for a nucleus were taken as a given quantity. There was no 
constructive description of how such structure functions arise. 

Recent work by McLerran and Venugopalan has given rise to a picture of how 
the structure functions arise at small x for very large nuclei at ultra-relativistic 
energy. In this description, the effects of quantum and charge coherence of the 
partons in the nuclear wavefunction are properly included. The gluons arise from 
the non-abelian Weizsacker-Williams fields generated by the color charges of the 
valence quarks. 

In this paper, we will extend the treatment of McLerran and Venugopalan from 
a description of a single nucleus to the collision of two nuclei. This work is in some 
sense an extension of early effort which were somewhat ad hoc to describe such 
collisions by classical fields |^, ^. We will see that in the region where most of 
the parton density sits, the gluon distribution function can initially be described 
by a classical field. These classical fields are to be interpreted as resulting from 
coherently superimposing large numbers of gluonic quanta. This way the classical 
description, wherever applicable will automatically incorporate coherence effects. 

The gluon field for a single nucleus arising in this way is a non-abelian Weiz¬ 
sacker-Williams field. At the initiation of the collision, the non-abelian Weizsacker- 
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Williams fields of the two nuclei play the role of boundary conditions for the time 
evolution of the gluon field. This classical field eventually evolves into gluon quanta. 

The picture we have of the collision is therefore the following. Before the nuclei 
collide, they are described by valence quarks and their coherent Weizsacker-Williams 
fields. These fields are classical in the sense of classical electromagnetic fields, but 
of course can not be thought of as composed of particle with classical phase space 
distributions. During the collision, the fields are still classical, but sufficiently strong 
so that the equations of motion evolve the fields non-linearly with time. As time 
evolves, the field weakens. When the strength of the gluon field is sufficiently low, 
the field equations linearize, and the gluon field describes the evolution of weakly 
interacting classical gluon waves. At this time, the coherent addition of the fields 
is no longer important, and they should be described by an incoherent distribution 
of gluons. The parton cascade model may therefore be used. 

Prior to this time however, the coherence in the gluon field is essential. The 
simple fact that the evolution of the gluons is described by a classical field is a con¬ 
sequence of the fact that the gluons are in some locally coherent state. A description 
in terms of incoherent classical particles is simply not possible. 

In the second section, we review the relevant results of computation of the 
small X structure functions for a single large nucleus. We will attempt to describe 
the kinematic limits of applicability of this description. We will argue that the 
Weizsacker-Williams fields should describe the distribution of gluons in the region 
of transverse momenta which gives the dominant contribution after integrating over 
transverse momenta. 

In the third section, we set up the problem of nucleus-nucleus scattering. We 
derive an equation for the time evolution of the gluon field. We relate the results of 
such a computation to the phase space density of gluon radiation. 

In the fourth section, we summarize our results and speculate on their region 
of validity. 
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2 Review of the McLerran-Venugopalan Model 


In the work of McLerran and Venugopalan Q, it was argued that for very large 
nuclei, ^ oo at small values of Bjorken x, x « the quark and gluon 

distribution functions are computable in a weak coupling limit. This is because the 
density of partons per unit area defines a dimensionful scale and when 


2 1 dN , 2 

-->> A-r 


( 2 ) 


7ri?2 dy ^ 

the strong coupling parameter as{p^) should become small. Here y ~ ln(l/x). 

In lowest order in a naive weak coupling expansion, it was shown that the gluon 
distribution function was of the Weizsacker-Williams form, that is proportional to 
1/x. It was also shown that the p±_ dependence was also of the Weizsacker-Williams 
form dN/d?p^ ~ ^/p\ for olsP « Pi. « P where p AqcdA^/^ 

Of course the naive weak coupling expansion may not be strictly valid, since 


there is the well known Lipatov enhancement of the low x structure functions |K]. 
This enhancement involves quantum corrections to the lowest order naive weak cou¬ 
pling result, and changes the small x distribution to While this behavior 

is computable in the McLerran-Venugopalan model, its nature is not yet fully under¬ 
stood. We expect however that as far as the local effects on the parton distribution 
at fixed rapidity, y ~ ln(l/x), the main effect is to renormalize the charge which 
generates the Weizsacker-Williams field. 

The charge which generates this field in lowest order in the naive weak cou¬ 
pling expansion is the charge of the valence quarks which are treated in a no recoil 
approximation. While it may be true that the Lipatov correction might involve 
new physics, and the picture might change, we will ignore its effects here except to 
state that we believe it will effectively renormalize the valence quark charge through 
some X dependent source of charge. To see how this might occur, recall that the 
strength of the Weizsacker-Williams distribution is proportional to the amount of 
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charge present at a value of x larger than that of the distribution. We are there¬ 
fore assuming the main effect of the quantum fluctuations is to generate an excess 
amount of charge at values of x larger than that at which we measure the parton. 

In the work which follow, we will not treat the problems generated by the Li¬ 
patov effect. We will instead concentrate on the naive lowest order approximation 
to the McLerran-Venugopalan model. This will be sufficient to understand many 
qualitative aspects of nucleus-nucleus collisions, and we hope in the end with small 
modifications can also be extended to include the effects of the Lipatov enhance¬ 
ment. 

In Refs. 1^, it was found that to compute the structure functions one simply 
treated the valence quarks as a source of light cone charge. 

Here the valence quarks are being treated as a source of charge moving at 
the speed of light along the light cone x~ = 0. The source of charge is being 
treated classically. This approximation as justified so long as the typical transverse 
momentum scale is 


pi_ « n (3) 

where p is proportional to the number of valence quarks per unit area. 

At the same time the number of gluon quanta at resolutions with p± « p 
will be sufficiently high to allow for a description of the gluonic degrees of freedom 
through a classical field. 

Within these limits, all one has to do is to formulate and solve the Yang-Mills 
equations in the presence of the classical current induced by the valence quarks: 

= U[A]ix,zix))rizix))U[A]{zix),x) . (4) 

Here z{x) = x\^+^q serves as a reference point used to define “initial values” for 
the color distribution of the valence quarks 

r{z{x)) = 6'^+6{x-)p{x^) (5) 
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which then, due to covariant current conservation [x)] = 0 evolve along 

the particles trajectory via parallel transport or link operators U\A]{x^ z{x)) := 
Pexp— A~{x'~^^x±,x~ = 0) connecting the points x and z{x) along the 
particles trajectory. 

Given a solution to the equations of motion, charge density p is to be treated 
as a stochastic variable, and to compute ground state expectation values one must 
average over all sources with a local Gaussian weight 

j[dp\ expj-y (fxx^p{xA_f^ (6) 

This Gaussian distribution arose from the approximations used in Ref. 0. It was 
argued there that on the transverse resolution scales corresponding to pT << P 
that the valence quark charges may be treated classically. The exponential factor is 
the contribution to the phase space density associated with counting the number of 
states of valence charges for a fixed value of the classical charge. It can be thought 
of as arising from the following classical picture. Suppose we look in a tube through 
the nucleus. This tube has a transverse size much less than a Fermi but large 
enough so that it intersects many nucleons. In this case, there will typically be 
many valence quarks inside the tube each coming from a different nucleon. The 
color charge of each quark is therefore uncorrelated with that of any other quark 
and the color charges will add together in a random walk. This will lead to the 
above Gaussian distribution. 

Physically, the picture one has is the following: The valence quarks are recoilless 
sources of color charge propagating along the light cone. Their charge can fluctuate 
from process to process and the averaging over charges corresponds to this fluctu¬ 
ation. The local charge density is therefore a random variable. The reason why 
such a stochastic source of charge arises is because the transverse resolution scales 
which we are interested in are small compared to a fermi. On such a scale, when one 
looks at the nucleus, one sees uncorrelated quarks coming from different nucleons. 


8 


The source of color charge therefore random walks in color space. The criteria that 
p± « is the criteria that within each transverse resolution scale, there are many 
quarks so that the color charge is typically large and can be treated classically. 

The solutions of the above Yang-Mills equations can be chosen to be of the form 

= 0 
A- = 0 

= e{x~)ai{x±) (7) 

Here the first line may be interpreted gauge choice. Using light cone gauge A'^ = 0 
one has direct access to the gluon distribution functions of the parton model. The 
requirement to have A~ = 0 then could still be implemented as a gauge choice at 
least along the trajectories of the particles, making use of the residual gauge freedom 
present in any axial gauge. In this case it turns out that there is a particular solution 
to the equations of motion which has A~ vanishing everywhere. On such a solution 
the link operators on the right hand side of the Yang-Mills equations drop out 
entirely and the equations become 

= 0 

V • a = p{x±) (8) 

The solution to these equations is that 

ai = --U{xx)yiU^ix±) (9) 

ig 

that is a pure two dimensional gauge transform of a the vacuum. 

Physically, this solution is also easy to understand. The solution is a gauge 
transform of vacuum on one side of the sheet of valence charge, and another gauge 
transform of vacuum on the other side of the sheet. We have chosen the field to be 
zero on one side of the sheet as an overall gauge choice. (This could be relaxed by 
an overall gauge transformation). Because of the discontinuity in the fields at the 


9 



sheet of valence charge, the solution is not a gauge transform of the vacuum fields. 
Its discontinuity gives the source of valence charge. 

Although we have not been successful in explicitly finding the solution to this 
equation, it is in principle possible to do numerically. Several generic features of 
the averaging over different sources of charge are possible to infer nevertheless. 

For p± < asp, the typical value of the external charge is so large that it is a 
bad approximation to linearize the gauge transformation and directly compute the 
held in a naive weak coupling expansion. In this region, the non-linearities of the 
held equation become important. 

In this kinematic region, the shape of the Weizsacker-Williams distribution 
changes form, as is shown in Fig. 2. 



Fig. 1: Weizsacker-Williams 
distribution for a single nucleus 


The l/p\ behavior turns over and goes to a constant. This provides a low 
momentum cutoff in the number of gluons generated by the distribution. At high 
momentum, we can compute no further than p± < p. The distribution should nev¬ 
ertheless extend beyond this region. In fact the upper momentum cutoff should be 
determined only by the kinematic limit of the process considered. Strictly speaking 
the number of Weizsacker-Williams gluons is inhnite, but only logarithmically, and 
the cutoff will be determined by the process of physical interest. 
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Suppose we consider the production of gluons in a two nucleus collision. We can 
naively determine the cutoff in momentum of produced gluons. Lets us assume that 
we are in the weak coupled perturbative region. The rate must be proportional to 
the density squared of gluons, so that it must be proportional to It involves 

scattering, so their are two more factors of as- Therefore on dimensional grounds, 
we expect that the probability of making a pair of on mass shell gluons should be 
of order /p\. This is of order one when p_i_ ~ ap, that is at precisely the place 
where the fields evolve non-linearly. It cuts off rapidly in p_L, so that the number of 
produced gluons should be ultraviolet finite. 

This example teaches us two things: First that for a physical process there is 
no ultraviolet divergence and for this process the important contribution for gluon 
production is at scales less than jj,. Second, that the process is strongest in the 
region where the field is strong. In this region, the field is evolving non-linearly, and 
the coherence of the field is important. It would therefore be a mistake to assume 
that the distribution of produced gluons reflects the distribution in the initial nuclei. 
This is true for gluons with p± > asp, that is ’hard gluons’, but the softer gluons 
which dominate the production are in a non-linear region. 

The appearance of these non-linearities might be qualitatively included in the 
parton cascade model. However insofar as there is an infrared cutoff dependence in 
some physical process, the results will be somewhat quantitatively unreliable. Pro¬ 
cesses without such a cutoff dependence would of course be more reliably computed. 

The hope will be in our attempt to compute gluon production is that the clas¬ 
sical non-linearities will cutoff the the naive divergence in the production amplitude 
at small p±. To see that this is plausible, recall that the single gluon distribution 
changes its from at small p± from 1/p^ to constant at p± < asp- It is therefore 
quite plausible that these effects in fact cutoff the singularity at some scale of order 
asp- If this is so, then if A is large enough so that >> 1, this cutoff is at a 
scale much larger than Aqcd, and the computation is self-consistent. 
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3 The Two Nucleus Problem 


We now turn to the problem of nucleus-nucleus scattering. We work in the center of 
mass frame. Both nuclei are at sufficiently high energy so that they can be treated as 
infinitesimally thin sheets. They are large enough so that these sheets can be taken 
to be of infinite extent in the transverse direction. We will be interested in describing 
the production of gluons at typical momentum scales which are p± « fJ-, but much 
larger than 200 MeV. In this case the source of color charge can be taken as classical 
as in the McLerran-Venugopalan model. The charge is of course a stochastic variable 
which must be integrated over with a Gaussian weight as described in the previous 
section. 



Fig. 2: Two nuclei Lorentz- 
contracted to infinitely thin 
sheets before the collision takes 
place 


The sources of color held set up a classical color held. After the collision the 
color held will begin to evolve in time. Much after the collision, the color held will 
describe the propagation of free gluons. In this section we will describe how to 
compute the evolution of the color held, and then how to compute the hnal state 
distribution of gluons. 

The Yang-Mills equation for the two source problem is 

( 10 ) 
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where 


J+ = 6{x~)pi{xi_) 

J- = 6{x'^)p2{Xi_) 

J* = 0 (11) 

and we have restricted ourselves to work in a gauge where the link operators along 
the particle trajectories drop out. 

Before the collision takes place, we find a solution of the equations of motion 
to be 


A+ = 0 
A- = 0 

= 0{x~)6{—x'^)a\{x±_) + 9{x'^)6{—x~)a\{x±) (12) 


This is a solution of the Yang-Mills equations in all of space-time except on or 
within the forward light cone, as shown in Fig. In the forward light cone, we 



Fig. 3: Regions with different 
structures of the gauge poten¬ 
tial: 

In regions 1 and 2 we have the 
well known one nucleus solu¬ 
tions ai^ 2 - While in the back¬ 
ward light cone there the gauge 
potential is vanishing we have 
a nontrivial solution in the for¬ 
ward lightcone, region 3 


must add in extra pieces in order to have a solution. 


This will be done below. The 


13 






two dimensional vector potentials are pure gauges and solve for t < 0 

V • ai,2 = Pi,2{xa_) (13) 

The physical picture one has of this solution prior to the collision is that the 
nuclei have zero held in front of them as they approach one another. Behind them 
the held is pure gauge. Because each nucleus has a different charge density, the 
gauge is different for each nucleus. This is an exact solution of the equations of 
motion so long as one is outside the forward light cone, that is in regions which are 
out of causal contact from the collision event. 

The fact that we have a solution of the equations of motion which does not 
evolve in time before the collision is remarkable. This solves the problem of cascade 
models that an isolated nucleus composed of partons will spontaneously fall apart. 
Here the individual nuclei and their parton clouds are static except for their overall 
center of mass motion. 

How do we describe the fields after the collision? Except at the forward light 
cone, that is, when we are inside or outside the cone, the fields satisfy free field 
equations. We will look for a Lorentz covariant solution to the equations of motion. 
We try, for > 0 and > 0 

= x^a(r, x_l) 

A~ = x“/3(r, x_l) 

= a\{T,XA_) (14) 

where 

r = 2x+x“ (15) 

This solution only depends on the longitudinal boost invariant variable r and has 
no dependence on the space-time rapidity variable 

r/=lln^ (16) 
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The factors of x"*" and x~ in the definition of the vector potential guarantee that 
under longitudinal boosts, the vector potential transforms properly. 

By making a gauge transformation which is only a function of proper time and 

x±, 


U = U{t,x±) 


(17) 


we see that we can fix 


P{r,x±) = -a{T,x±) 


(18) 


which we shall choose to do. This choice corresponds to a gauge condition 


x^A + x A'^ = 0 


(19) 


which in turn is consistent with dropping the link-operators in (|T^). 

If such a solution solves the equations of motion an boundary conditions, it will 
predict that the distribution of partons is boost invariant. It is the generalization 
therefore of Bjorken’s boost invariant hydrodynamic equations to the equations 
which generate the initial conditions for the hydrodynamic equations. 

As a consequence of the boost invariance of the Yang-Mills equations, the ansatz 
above solves the equations within the forward light cone. This can be checked 
explicitly. The equations which result for 0 : 3 , and a, for x^,x~ >0 are 


■^drT^dra — 
T'* 

1 
T 




D\a 


1 


drTdrol — igr"^ 


Di, dra\ + igr [a, dr a] 
a, \D\a\] - 


0 

0 

0 


( 20 ) 


These four equations can be checked to be consistent with on another. Prom 
this point on, all vector indices will refer to two dimensional transverse vectors. The 
longitudinal and time coordinates will be denoted separately. 
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In intermediate step in deriving the above equations, we computed the field 
strengths . The results of that computations are 

F~^~ = ——drT^a 

T 

ai — a\ — ig a\,ai 

= -x^ (^^drai^[D\a\^ ( 21 ) 

The only task which remains is to show that the above solution also satishes 
the boundary conditions generated by the sources. For either < 0 or x"*" < 0, 

the ansatz above satishes the equations trivially. We look hrst at the equation 
[D^, F^*] = 0. This has a delta function singularity at x“ = x"*" = 0 which requires 
that 

al{0,x±) = a\{x±) + al{x±) ( 22 ) 

r=0 

and there are no further discontinuities in this equation. 

Now for [D^, F^^] = J^, we hnd that 

^9 r ' 1 

a{T,x±) ^ = -17 aii{x±),a2ix±) ( 23 ) 

T=0 z L -I 

These two equations for reduce to the same boundary condition, and therefore 
neither a nor a\ are over constrained, demonstrating once more that our ansatz 
contains the correct degrees of freedom. 

Note that assuming the boundary conditions above, we are implicitly requiring 
that the solution be regular at r = 0. It is easy to check that the quantities a 
and can either be regular at the origin or diverge like a ~ 1/r^ and ct ~ ln(r). 
These singular solutions will lead to a divergent energy density, and are therefore 
not allowed. 

So with the above two boundary conditions, the solution to the equations of 
motion are uniquely specified. This solution is remarkable since in spite of the 
possible asymmetry in the charge on either nuclei, the solution is up to trivial 
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factors rapidity independent. This has amusing phenomenological consequences for 
the collisions of asymmetric nuclei. The distribution should be flat in rapidity. The 
height of the central plateau of course depends on the asymmetry in a non-trivial 
way. 

In order to determine the gluon radiation produced by these fields, we must 
solve them at proper times long after the collision. We expect that the energy 
density will dissipate and therefore the field strengths will become small. Using the 
expressions above, we conclude that asymptotically, for large r, 

a{T,x±) V{xi_y{T,Xi_)V\xx_) 

Oil{T,Xi_) V{x±) €\t,x^) - v\x±) (24) 

The solution should tend to a small field plus a gauge transformation. The value of 
this gauge transformation is determined by the field equations and has a non-trivial 
dependence on the sources. It results from solving the non-linear time evolution 
equations for the fields. 

The equations of motion for the fields in the asymptotic region are linear for e 
and e*. The equations are 

J-drT^dre-V^e = 0 

^drTdre^ - e^' = 0 (25) 

Observe, that V*e* does not enter the asymptotic equations, so that there are in 
fact only two dynamical degrees of freedom in the solution, as must be the case. 
The solutions to the above equations at asymptotically large r are of the form 

In this equation, the frequency u; =| kx_ |, and the vector 

K^ = e^^yiuj (27) 
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The notation +C.C. means to add in the complex conjugate piece. 

To derive an expression for the energy density, we recall that r is large. Near 
z = 0, this implies that the range of z where t ^ t » z that the solutions are 2 
independent. This means they asymptotically have zero pz- Now suppose we are at 
any value of z, and r is large but t ^ z. We can do a longitudinal boost to z = 0 
without changing the solution. Again in this frame the solution has zero pz- We 
see therefore that for the asymptotic solutions that the space time rapidity is one 
to one correlated with the momentum space rapidity, that is at asymptotic times 
we find that 


r/=iln(x+/x )=y=^ln{p^/p ) 


(28) 


To proceed further, we compute the energy density in the neighborhood of 
z = 0. Here asymptotically t = t. The energy in a box of size R in the transverse 
direction and dz in the longitudinal direction, with L « t becomes |13] 


dE = 


dz f d^k\ 


t 


(2vr) 


tLO 


T. I 


(28) 


2,6 


Recalling that dy = dz/t, we find that 

dE 1 


dyd?k_\_ (27r)3 


OJ 




(30) 


2,6 


and the multiplicity distribution of gluons is 


dN _ 1 dE 
dyd‘^k± uj dyd?k± 

As we expect for a boost covariant solution, the multiplicity distribution is rapidity 
invariant. 

Finally, we must comment a bit on the characteristic time scale for the dissi¬ 
pation of the non-linearities in the equations for the time dependent Weizsiicker- 
Williams fields. This is difficult to estimate in general, but scaling arguments should 
suffice to estimate the time scale. The basic point is that the typical momentum 
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scale in the problem relevant for the formation of most of the gluons is p± ~ ap 
which up to logarithms is p± ~ 20QMeV. The characteristic time scale for the 
dissipation of the classical non-linearities should therefore be of order r 1/P±- 
This is in agreement with other estimate of the characteristic formation time for 
partons, and is in agreement with the model of Geiger and Muller. 

4 Summary and Conclusions 

We have derived a theory of the formation of gluons which is applicable for small 
X gluons in the collisions of very large nuclei. We have found that the gluon dis¬ 
tributions as measured in deep inelastic scattering undergo an entirely non-trivial 
evolution in forming gluons which would be the initial conditions for a parton cas¬ 
cade. We have shown how one can compute these initial conditions. 

After finding the initial conditions, the subsequent evolution might be described 
by a combination of the parton cascade and hydrodynamics. 

There are many further problems to be addressed in this theory. The equations 
described above must be numerically solved. This will provide for initial conditions 
for average head on collisions. In addition it will predict the spectrum of fluctuations 
from collisions to collision. Perhaps the most interesting problem is to compute the 
hard particles produced during the early evolution of the distributions so as to find 
a precise quantitative test of the theory. 
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